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\0 ' Abstract 



The formalism of photon and dilepton productions in a ther- 
mal bath are discussed in the frame work of finite temperature 
field theory. The emission rate has been expressed in terms of the 
Tlj" . discontinuities or imaginary part of the photon self energy in the 

I thermal bath. The photon and dilepton emissions have been com- 

' puted by taking into account the in-medium modifications of the 

On . particles appearing in the internal loop of the self energy diagram 

using a phenomenological effective lagrangian approach. The re- 
sults has been compared with the photon spectra measured by 
WA 98 collaboration at CERN SPS energies. 



'. I. Introduction 

> ■ 

•i-H . 

^ ' Numerical simulations of the QCD (Quantum Chromodynamics) equa- 

■ tion of state on the lattice predict that at very high density and/or tem- 

perature hadronic matter undergoes a phase transition to Quark Gluon 
Plasma (QGP) [|ll(see Fig. p. One expects that ultrarelativistic heavy 
ion collisions might create conditions conducive for the formation and 
study of QGP. Various model calculations have been performed to look 
for observable signatures of this state of matter. However, among various 
signatures of QGP, photons and dileptons are known to be advantageous, 
primarily so because, electromagnetic interaction could lead to detectable 
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Early Universe 




Baryon density p p/p„ = 5-10 



Figure 1: QCD phase diagram. 

signal. However, it is weak enough to let the produced particles (real pho- 
tons and dileptons) escape the system without further interaction and 
thus carrying the information of the constituents and their momentum 
distribution in the thermal bath. 

The disadvantage with photons is the substantial background from 
various processes (thermal and non-thermal) 0. Among these, the con- 
tribution from hard QCD processes is well understood in the framework 
of perturbative QCD and the yield from hadronic decays e. g. 7r° — > 7 7 
can be accounted for by invariant mass analysis. However, photons from 
the thermalised hadronic gas pose a more difficult task to disentangle. 
Therefore it is very important to estimate photons from hot and dense 
hadronic gas along with the possible modification of the hadronic prop- 
erties. 

The importance of the electromagnetic probes for the thermodynamic 
state of the evolving matter was first proposed by Feinberg in 1976 
While for most purposes one can calculate the emission rates in a classi- 
cal framework, Feinberg showed that the emission rates can be related to 
the electromagnetic current current correlation functions in a thermalised 
system in a quantum picture and, more importantly, in a nonperturba- 
tive manner. Generally the production of a particle which interact weakly 
with the constituents of the thermal bath ( the constituents may interact 
strongly among themselves, we need not give any explicit form of their 
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coupling strength) can always be expressed in terms of the discontinu- 
ities or imaginary parts of the self energies of that particle Q which 
never thermalizes. We should, therefore look at the connection between 
the electromagnetic emission rates (real photons and lepton pairs, which 
correspond to virtual photons) and the photon spectral function ( which 
is connected with the discontinuities in self energies) in a thermal sys- 
tem 1^, which in turn connected to the hadronic electromagnetic current 
current correlation tensor [0 through Maxwell equations. We will show 
below that the photon emission rate can be obtained from the dilepton 
(which originate from a virtual photon) emission rate with little effort. 

In section II we present the general formalism for the dilepton and 
photon production rate from a thermal bath. In section III we calculate 
the medium modifications of hadrons. Section IV is devoted to discuss 
the results of our calculations. In section V we present a summary and 
discussion. 

II. Electromagnetic Probes - Formulation 

We begin our discussion with the dilepton production rate. In the follow- 
ing we will follow the work of Weldon for the derivation of dilepton 
emission rate. Let Cq Jfj, be the interaction between the photon and 
the particles in the heat bath, then to lowest order in the electromag- 
netic coupling, the S matrix element for the transition | /) — > | if; 
is given by 

SHi = eo J {H-l+r I Jf.A^'ix) I /)rf^xe^^-^ (1) 

where | /) is the initial state corresponding to the two incoming nuclei, | 
H] is the final state which corresponds to a lepton pair plus anything 
(Hadronic), the parameter cq is the bare (unrenormalised) charge and 
q = Pi + P2 is the four momentum of the lepton pair. Since we assume 
that the lepton pair do not interact with the emitting system, the matrix 
element can be factorised in the following way, 

{H-l+l- I J^A^'ix) \ I) = {H\ A'^ix) I /)(/+/- I J,, I 0) (2) 
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Putting the explicit form of the current in terms of the Dirac spinors 
we obtain 

S., = e.|^^/d'.e-(i/M^W|/> (3) 

Therefore dilepton multiphcity from a thermal system is obtained as, (by 
summing over the final states and averaging over the initial states with 
a weight factor Z{I3)-^ e'^^') 

^—EElS^.re-' (4) 

where Ej is the total energy in the initial state. After some algebra N 
can be written in a compact form as follows, 

where L^,^ is the leptonic tensor given by 

spins 

= P^P2+P^P2-y^'^^ (6) 
and Hj^i, is the photon tensor 

(7) 

to obtain the above equation we have used the resolution of identity 
1 = J2i I I) {I I ; and the energy conservation equation Ej — Eh + qo 
where Qq is the energy of the lepton pair and Eh is the energy of the rest 
of the system produced after collision. Using the translational invariance 
of the matrix element we can write 

//^, = 27rQe-^>'^'^ (8) 

where is the four volume of the system and p'^'^ is the spectral function 
of the photon in the heat bath, 

P,M<lo) = ^^r^ / d'xe'^-Y.{H I A,{x)A^{0) \ H)e-^^- (9) 



H 
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The rate of dilepton production per unit volume (N/Q) is given by, 



— - zvr CqJ.^^ e 



{2'Kf2Ei (27r)32E2 



(10) 



This result which expresses the dilepton emission rate in terms of the 
spectral function of the photon in the medium is a fundamental result. 
At zero temperature {(3 —>■ oo) the only state which enters in the spectral 
function (p(s)) is the vacuum and the spectral function can be expressed 
in terms of the imaginary part of Greens function (i.e the discontinuity 
along the real axis) p{s) = Im G'(s) p. 

Unlike at zero temperature, in the case of non-zero temperature the 
time ordered {D^") and the Feynman propagator {Dp') are different. 
The relation between them is given by H], 



Im D 



:i + 2/BE)ImL> 



(11) 



where Jbe is the thermal distribution for photons. Now the time ordered 
Green's function can be expressed in terms of the spectral function as 
follows 



da 



le 



— a + It 

Using the KMS relation one can show that 

substituting Eq. (p!3| ) in Eq. (|T2|) and using Eq. ( [Tl| ) we obtain 

-1 



e^'^p'^i-q^.-q) 



P^'iqo^q) 



TT 



;i + /B£;(go))Im£T(9o, q) 



(12) 



(13) 



(14) 



The above equation implies that to evaluate the spectral function at 
T 7^ we need to know the imaginary part of the Feynman propagator. 
It is important to note that above expression for spectral function reduces 
its vacuum value as /3 ^ oo. The Feynman propagator can be expressed 
in terms of the proper self energy as 



g2 -Ht g2 _ n, 



+ (C-i 



q^iqu 



(15) 
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where is the wave function renormahzation constant at zero temper- 
ature, A^^ and B^" are transverse and longitudinal projection tensors 
respectively. The presence of the parameter C, indicates the gauge de- 
pendence of the propagator. Although the gauge dependence cancels out 
in the calculation of physical quantities, one should, however, be careful 
when extracting physical quantities from the propagator directly, espe- 
cially in the non-abelian gauge theory. 
Using Eqs. (0), (|l|) and (|15D we get 

§- . 2e^L,M-r. + (16) 
where = Z^e^, is the physical charge of the electron and 

" {q^ -ReUy + lmUl ^ 

The emission rate of dilepton can also be expressed in terms of the 
electromagnetic current-current correlation functions 0. Denoting the 
hadronic part of the electromagnetic current operator by J^, the leptonic 
part by Jl and the free photon propagator by D'^'^, we have the matrix 
element for this transition : 

Shi = {H-J+l- \J d''xdSj''^{x)D^' JUv) I /) (18) 

As in the earlier case the leptonic part of the current can be easily factored 
out. Writing the Fourier transform of photon propagator and squaring 
the matrix elements, one obtains, for the rate of dilepton production, 

where W^'^{q) is just the Fourier transform of the thermal expectation 
value of the real time electromagnetic current- current correlation function 

/p—PEjj 
d^xe-"^''J2(H I J''{x)r{0) I H)^— (20) 
H ^ 

The subtleties of the thermal averaging have been elucidated earlier. It 
is thus readily seen (eq. |19D that the dilepton data yields considerable 
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information about the thermal state of the hadronic system; at least 
the full tensor structure of W^'^ can in principle be determined. The 
most important point to realize here is that the analysis is essentially 
nonperturbative. 

At this point one may wonder what is the connection between the 
electromagnetic current-current correlation function and the spectral func- 
tion? The connection can be expressed in a straight forward way by sub- 
stituting and J,y using Maxwell equation {dad" — 9^' {da A"") = j^) 
in Eq. ( pO|) to obtain, 

el = 27r (g2 g^^- - (^2 _ (21) 

In the literature most of the dilepton production rate from a thermal 
system are calculated with the approximation 11 j- = I\.l which gives the 
spectral function as 



P^''{qo,q) = {-g^' + q^q'/q')—{l + fBE{qo))l^ 

71 



1 



q'-T^I 



(22) 



Since Imllfc and Rellfc are proportional to a (the fine structure constant) 
they are very small for most of our practical purposes. Therefore pk is 
given by 

Pk = — — 23 
Using Eqs. (0), (0), (^) and the following result 



^C^- (24) 
q2 



where C^*^ = q^^ — q^ g'^'^, we get, 



dR a , 2m , / 4m^ „ , , , , 



This is the familiar results most widely used for dilepton emission rate [|]. 

To obtain real photon emission rate per unit volume {dR) from a sys- 
tem under thermal equilibrium we note that the dilepton emission rate 
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differs from tlie plioton emission rate in the following way. The factor 
L^y is the product of the electromagnetic vertex 7* l~ , the 

leptonic current involving Dirac spinors and the square of the photon 
propagator should be replaced by the factor J2 ^^i^vi= ~g^iv) for the real 
(on-shell) photon. Finally the phase space factor d?'pi/[{2TxY2 E]\ (i^p2/[(2vr)^2£'2 
should be replaced by, (i'^g/[(27r)^2go] to obtain 

The current-current correlation function is related to the photon self 
energy as |l^, [ill] 



W^M- = 2/Bi5(go)Imn^, (27) 

Therefore, 



The above equation can also be obtained directly from eq. The 
emission rate given above is correct up to order in electromagnetic 
interaction but exact, in principle, to all order in strong interaction. But 
for all practical purposes one is able to evaluate up to a finite order of loop 
expansion. Now it is clear from the above results that to evaluate photon 
and dilepton emission rate from a thermal system we need to evaluate the 
imaginary part of the photon self energy. The Cutkosky rules at finite 
temperature [0, 0, [l^ gives a systematic procedure to calculate the 
imaginary part of a Feynman diagram. The Cutkosky rule expresses the 
imaginary part of the n-loop amplitude in terms of physical amplitude 
of lower order {n — 1 loop or lower). This is shown schematically in 
Fig. ^. When the imaginary part of the self energy is calculated up to 
and including L order loops where L satisfies x + y < L + 1 then we 
obtain the photon emission rate for the reaction; x particles — > y particles 
+ 7 and the above formalism becomes equivalent to the relativistic kinetic 



theory formalism [p!5|]. 



For a reaction 1 + 2 8+7 the photon (of energy E) emission 
rate is given by |TB[, 



E— = — — — / ds / dt M ' / dEi 

d^p 16(27r)7E i(mi+m2)2 ' I I / 1 
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Figure 2: Optical Theorem in Quantum Field Theory 
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In a similar way the dilepton emission rate for a reaction a a 
can be obtained as 



dN 



d'^xd^qdM 



J 2£;„(27r)3-^^^"^ 



d^Pa 



2E-a{2ny 



-J'iPa) 



l+l- 
d'p2 

2Ei(27r)^J 2£;2(27r) 



d'pi 
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I M {2T^)%Pa +p-a-pi- P2)5{Ea + E, - E, - E^) 

6{pa + p-a)S{M-E,-E-,) (30) 

where = {pa +Pa)^ is the invariant mass and f{pa) is the occupation 
probabihty in the momentum space. The Pauh blocking of the lepton 
pair in the final state has been neglected in the above equation. 

III. Finite Temperature Properties 

To calculate the effective mass and decay widths of vector mesons we 
begin with the following VNN (Vector - Nucleon - Nucleon) Lagrangian 
density: 

CvNN = gvNN [n^^t'^NV,^ - ^Na,,T^Nd''V^^ , (31) 

where = {uj^, M is the free nucleon mass, is the nucleon field 
and Ta = {1, t}. The values of the coupling constants gvNN and will 
be specified later. With the above Lagrangian we proceed to calculate 
the p-self energy. 

d'^p 



/dp 
j^K^^{p,k), (32) 

with, 

Tr [r,{k) {i> + M*)r,(-fc) {i>-lf: + M*)] 
(p2_M*2)((p-A;)2-M*2) ' 

The vertex T^{k) is calculated by using the Lagrangian of Eq. (^) and 
is given by 

r,(A:) = 7. + (34) 

Here M* is the in-medium (effective) mass of the nucleon at finite tem- 
perature which we calculate using the Mean-Field Theory (MFT) [p!7[] . 
The value of M* can be found by solving the following self consistent 
equation: 

+ (35) 
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where /Ar(T)(/^(T)) is the Fermi-Dirac distribution for the nucleon (antin- 
ucleon), rris is mass of the neutral scalar meson (a) field, and, the nucleon 
interacts via the exchange of isoscalar meson with coupling constant Qs- 
Since the exact solution of the field equations in QHD is untenable, these 
are solved in mean field approximation. In a mean field approximation 
one replaces the field operators by their ground state expectation val- 
ues which are classical quantities; this renders the field equations exactly 
solvable. 

The vacuum part of the p self energy arises due to its interaction 
with the nucleons in the Dirac sea. This is calculated using dimensional 
regularization scheme (see |jT8| for details). 

In a hot system of particles, there is a thermal distribution of real 
particles (on shell) which can participate in the absorption and emission 
process in addition to the exchange of virtual particles. The interaction of 
the rho with the onshell nucleons, present in the thermal bath contributes 
to the medium dependent part of the p-self energy. This is calculated 
from Eq.(p2]) using imaginary time formalism as. 



Reflmprifcj, k 



0) 



dp 



X 



+ 



^(2p2 + 3M*2)+^2|^*.'^v^ 



(36) 



where, uj"^ = p^ + M 



*2 



IV. Results 

In this section we present the results of our calculation of the effective 
masses and decay widths of vector mesons and its effect on electromag- 
netic probes. For our calculations of p and uj mesons effective masses we 
have used the following values of the coupling constants and masses : 
= 6.1, glNN/4:7i = 0.55, m,= 550 MeV, nip = 770 MeV, M = 938 
MeV, g'^jAn = 9.3, glNN/4:TT = 20, and = 0. In Fig. 3 we show the 
variation of the ratio of effective mass to free mass for hadrons as a func- 
tion of temperature. The expected trend based on Brown - Rho scaling 
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is also shown. Firstly, we do not observe any global scaling behaviour. 
Secondly, in our case the effective mass as a function of temperature falls 
at a slower rate. 

Many authors have investigated the issue of temperature dependence 
of hadronic masses within different models over the past several years. 
Hatsuda and collaborators ||2y, ^ and Brown showed that the 



use of QCD sum rules at finite temperature results in a temperature 
dependence of the qq condensate culminating in the following behaviour 
of the p-mass: 

^-(-^)" (37) 

nip V T^J 

where Ty is the critical temperature for chiral phase transition. Brown 



and Rho |]2^ also showed that the requirement of chiral symmetry yields 



an approximate scaling relation between various effective hadron masses, 

(38) 



M nip 

These calculations show a dropping of hadronic mass with tempera- 
ture. Calculations with non-linear a-model, however, predict the op- 



posite trend |^| 
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Figure 4: Rho decay width (Fp) as a function of temperature (T). 
Dashed and sohd hues show calculations of rho width with effective mass 
due to nucleon loop, with and without BE. Dot dashed line represents 
the same but with effective mass due to pion loop. 




The variation of the in-medium decay width (Fp) of the rho meson 
with temperature is shown in Fig. @. As discussed earlier, the effective 
mass of the rho decreases as a result of NN polarisation. This reduces 
the phase space available for the rho. Hence, we observe a rapid decrease 
in the rho meson width with temperature (solid line). However, the 
presence of pions in the medium would cause an enhancement of the 
decay width through induced emission. Thus when Bose Enhancement 
(BE) of the pions is taken into account the rho decay width is seen to 
fall less rapidly (dashed line); such a behaviour is observed quite clearly 
26[] . For the sake of completeness, we also show the variation of rho 



m 



width in the case where the rho mass changes due to nn loop. In this 
case, since the rho mass increases (though only marginally), the width 
increases (dot-dashed line). 

Now we present our results on photon emission rates from a hot 
hadronic gas. As discussed earlier, the variation of hadronic decay widths 
and masses will affect the photon spectra. The relevant reactions of 
photon production are tttt p'~f,'n-p ^7r7,7r7r —>■ 777, andvrr/ ~^ rj'j 
with all possible isospin combinations. Since the lifetimes of the rho and 
omega mesons are comparable to the strong interaction time scales, the 
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Figure 5: Total photon emission rate from hot hadronic matter as a 
function of photon energy at T=160, 180 MeV. The sohd and dot-dashed 
hues show results with and without in- medium effects respectively. Inset: 
Total photon emission rate is plotted in linear scale as a function of 
photon energy in the kinematic window, = 2.5 to 3.0 GeV at T =180 
MeV. 



decays p ^ vr tt 7 and a; — >• 7r° 7 are also included. The effect of finite 
resonance width of the rho meson in the photon production cross-sections 
has been taken into account through the propagator. 

In Fig. d^) we display the total rate of emission of photons from hot 
hadronic gas including all hadronic reactions and decays of vector mesons. 
At T =180 MeV, the photon emission rate with finite temperature effects 
is a factor of ~ 3 higher than the rate calculated without medium effects. 
At T = 160 MeV the medium effects are small compared to the previous 
case. 

The transverse momentum distribution of photon has been obtained 
by convoluting the production rate with the space time dynamics of the 
system. The boost invariant hydrodynamical model of Bjorken has 
been used for the space time evolution of the system with appropriate 
modification due to the shift in hadronic masses [^. In Fig. ^ we com- 



pare our results of transverse momentum distribution of photons with the 



preliminary results of WA98 Collaboration |29]. The experimental data 



represents the photon spectra from Pb -|- Pb collisions at 158 GeV per 
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Figure 6: Total thermal photon yield in Pb + Pb central collisions at 158 
GeV per nucleon at CERN SPS. The long-dash line shows the results 
when the system is formed in the QGP phase with initial temperature 
Tj = 180 MeV at Tj = 1 fm/c. The critical temperature for phase tran- 
sition is taken as 160 MeV. The solid (short-dash) line indicates photon 
spectra when hadronic matter formed in the initial state at Tj = 230 
MeV (Tj = 270 MeV) at Tj = 1 fm/c with (without) medium effects on 
hadronic masses and decay widths. 

nucleon at CERN SPS energies. The transverse momentum distribution 
of photons originating from the 'hadronic scenario' (matter formed in the 
hadronic phase) with (solid line) and without (short-dash line) medium 
modifications of vector mesons outshine the photons from the 'QGP sce- 
nario' (matter formed in the QGP phase, indicated by long-dashed line) 
for the entire range of pt- Although photons from 'hadronic scenario' 
with medium effects on vector mesons shine less bright than those from 
without medium effects for > 2 GeV, it is not possible to distinguish 
clearly between one or the other on the basis of the experimental data. 

So far we have studied the effects of in-medium hadronic masses and 
decay widths on the photon spectra, where, the latter is found to have 
negligible effects. However, the modification in the rho decay width due 
to BE, arising from the induced emission of pions in the thermal medium, 
plays a crucial role in low mass lepton pair production. It has been shown 
in Ref. that this effect arises naturally in a calculation based on finite 
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Figure 7: Dilepton emission rates as a function of invariant mass at T 
= 180 MeV and twice nuclear matter density. Solid and dotdashed lines 
show the results with and without BE effects respectively. The dotted 
line shows the result when no medium effect is considered. 

temperature field theory. 

For the sake of illustration, we consider the effect of thermal nucleon 
loop on the lepton pair production from pion annihilation via (rho me- 
diated) vector meson dominance. In the work of Li, Ko and Brown [pl[1 , 
the observed enhancement on the low mass dilepton yield was attributed 
solely to the the matter induced mass modification of the meson, while 
neglecting the BE effect. In Fig. (|^) the dilepton emission rate is plotted 
as a function of invariant mass at a temperature T=180 MeV. It is clear 
from the figure that the inclusion of BE effects reduces the yield by in- 
creasing the decay width. Quantitatively, a suppression by a factor of 3 
is observed at the invariant mass, M = m*. 

V. Summary and Discussions 

We have discuss the general formalism to evaluate the emission of elec- 
tromagnetic ejectiles from a thermal bath. This formalism has been 
applied to evaluate the photon and dilepton emission rate from a hot 
and dense hadronic system. We have calculated the effective masses and 
decay widths of vector mesons propagating in a hot medium. We have 
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seen that the mass of rho meson decreases substantially due to its in- 
teractions with nucleonic excitations and it increases only marginally (~ 
10-15 MeV) due to p — tt interactions. The overall decrease in the effec- 
tive mass is due to fluctuations in the Dirac sea of nucleons with mass 
M* and the in-medium contribution increases very little (~ 5 - 10 McV) 
from its free space vahie. The omega mass drops at T =180 MeV by 
about 250 MeV from its free space value. 

We have evaluated the rate of photon emission from a hadronic gas 
of TT, p, uj and r] mesons. It is seen that the photon rate increases by a 
factor ~ 3 due to the inclusion of in-medium masses at T =180 MeV. We 
observe that the inclusion of in-medium decay widths in vector meson 
propagator has negligible effects on photon emission rates, although the 
effect of in-medium decay widths with BE is substantial for the dilepton 
yield. 

We have compared the experimental data on photon spectra from 
Pb + Pb collisions at energies 158 GeV per nucleon with different initial 
conditions. We observe that photons from hadronic scenario dominates 
over the photons from QGP scenario for the entire pt domain. But in 
the hadronic scenario the photon spectra evaluated with and without in- 
medium properties of vector mesons describe these data reasonably well. 
Hence the transverse photon spectra at present do not allow us to decide 
between an in-medium dropping mass and a free mass scenario. Consid- 
ering the experimental uncertainty, it is not possible to state, which one, 
between the two is compatible with the data. Experimental data with 
better statistics could possibly distinguish among various scenarios. 
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